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ABSTRACT 

The evolutionary path followed by a galaxy shapes its internal structure, and, in particular, its 
internal colour variation. We present a study of the internal colour variation within galaxies 
from the Seventh Data Release of the Sloan Digital Sky Survey (SDSS DR7). We statistically 
study the connection between the internal colour variation and global galactic properties, look- 
ing for hints of the recent galactic evolution. Considering only galaxies with good photometry 
and spectral measurements, we define four luminosity-threshold samples within the redshift 
range 0.01 < z < 0.17, each containing more than 48000 galaxies. Colour gradients are cal- 
culated for these galaxies from the surface brightness measurements provided by the SDSS 
DR7. Possible systematic effects in their determination have been analysed. We find that, on 
average, galaxies have redder cores than their external parts. We also find that it is more likely 
to find steep colour gradients among late-type galaxies. This result holds for a range of classi- 
fications based on both morphological and spectral characteristics. In fact, our results relate, 
on average, steep colour gradients to a higher presence of young stars within a galaxy. Our 
results also suggest that nuclear activity is a marginal driver for creating steep colour gra- 
dients in massive galaxies. We have selected pairs of interacting galaxies, with a separation 
of 5 arcsec, in projected radius, and a difference in redshift of lOOkm/s, finding that they 
present steeper gradients than the average population, skewed towards bluer cores. Our analy- 
sis implies that colour gradients can be useful for selecting galaxies that have suffered a recent 
(minor) burst of star formation. 

Key words: galaxies: fundamental parameters, statistics, structure 



1 INTRODUCTION 



The subtleties of galaxy formation remain to be understood. The 
study of the internal structure of galaxies can help us to under- 
stand how they formed and evolved. Here we study the internal 
colour variation across galaxies, quantified through their colour 
gradient. Interest in internal colour variation has grown enormously 
in the last five years, especially due to the availability of both shal- 
low large survey s such as the Sloan Digital Sky Survey (SDSS, 
lYork et alJ bOQO). that can provide statistically significant samples 
of galaxies, and deep surveys, that provide high resolution images 
of galaxies at z > 0.5, suc h as those performed w i th the H ubble 
Space Telescope (HST; e.g. iMenanteau etHl]|2005l ; iFerreras et al.l 
120091) . Previous work has used the colour gradient as a measure of 
the galactic internal colour variation, which has proved to be very 
useful for exploring galactic properties, due to its tight dependency 
on the galactic distribution of metallicity, stellar age and dust (see 
e.g.lWu et al.l2005l ; lMicharj2005l ; lRawle et al.l201(il : lTortora etail 
2010). This, together with the accessibility of gradient measure- 
ments, since only photometry is needed, makes the colour gradient 
a useful tool to obtain statistical results for wide samples of galax- 



ies, that will be interesting to relate to different theoretical aspects 
from models of galaxy formation. 

In general, galaxies t end to have cores t hat are redder than 
their external parts (e.g. IPeletier et al.l Il990h . In this work we 
study the variation within nearby galaxies of the following opti- 
cal colours: (g - i), (g - r), (r - i), (i - z) and (r — z). Both (g - r) and 
(g - ;') rest frame colours are most sensitive to the intensity of the 
spectral break that appears around 4000 A, due to the accumulation 
of resonant absorption lines from metals in different stages of ion- 
isation (e.g. lHamiltor]|l985h . This break correlates with the tem- 
perature of the dominant stellar population, increasing for cooler 
stars, but it does not depend strongly on the galactic metallicity. 
Unlike late-type galaxies, early-type ones are dominated by K-type 
stars, showing a rather steep change in flux from g to r or i-band. If 
we consider that these colours are mainly related to the age of the 
galaxy then the redder core will imply that it is also older. 

The other three colours, (r - i), (r - z) and (i - z), in the rest 
frame directly depend on the proportion of main se quence to red gi- 
ant branch or older stars (Bruzual & Chariot 1993'). Thus, we could 
simply record these colours as depending on the relative amount of 
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old and young stars in the galaxies, though this is less clear for (i— z) 
and has a smaller dynamic range in the case of the (r - i) colour. 
Thus, we observe the same tendency as before: redder cores imply 
that in general galaxies tend to have older stars in their centres. 

However, these implications should be taken with caution for 
two main reasons: i ) Global colours can be affected by marginal 
young populations aLi & Had 120070 . ii) Colours within different 
galactic annuli do not necessarily agree well with global ones. Be- 
sides, broad band colours show a degeneracy in t heir dependence 
with stellar population age and metallicity (e.g. IWorthevI 1 19941 : 
Chario t. Worthev & Bressan 199rj j). Indeed, several studies (see e.g. 
Tamura & Ohtal |2004 IWu et all l2005l : lLa Barbera & de Carvalhol 
2009; Tortora et alj|2010l) have shown that metallicity and not age 



gradients are the principal origin of the colour variation. In the 
case of early-type galaxies, age gradients have been measured to 
be almost flat or positive while metallicity gradients are nega- 
tive, i.e. early-type g alaxies tend to have younger , higher metal- 
licity cores (see e.g.lSanchez-Blazquez eta l. 2006; Cleme ns" et al.l 
120091 : lLa Barbera & de Carvalhd 120091 : iRawle et alj|20ld) . Late^ 
types galaxies are more complex, showing a broad variety of age 
and also metallicity gradients, due to the interplay between dust 
content, stellar feedback, bursts of formation and stellar migration 
due to the inst abilities that trigger the formation of bars but also 
the spiral arms. iMacArthur etall d2009h observed that spiral galax- 
ies have bulges with positive age gradients, while their disks show 
negative ones, but there are examples of spiral galaxies with nega- 
tive age gradi ents for most of th eir extension, except for the most 
external radii dVlaiic et alj2009l) . 

Metallicity gradients within galaxies can be produced 
through d ifferent mech anisms. Winds generated by su- 
pernovae dLarsor] 1974) or active galactic nuclei (AGN) 
(Begel man. de Kool and Sikord [l991) can provide the metal- 
licity enrichment of the inter stellar or even the inter galactic 
medium. Due to the multiple phases that compose this medium 
the enrichment does not necessary happen in a symmetrical way, 
which can produce strong metallicity gradients and thus colour 
gradients. Galaxies with a significant age gradie nt can present 
very shallow or even inverted colour gradients dHinklev & Iml 
l200lh . since age gradients usually arise from localised regions of 
star formation, generally resulting in galaxies with bluer cores 
than their external parts. The presence of a diffuse component of 
dust could also explain some colour gradients. However, this has 
a marginal impac t, specially for early-type galax ies, which are 
almost dust free dMichard 120051 ; IWu et al.l l2005h . In summary, 
colour gradients in galaxies are likely to be related to the different 
star formation histories followed, which depend on the way gas is 
distributed and collapses after accretion and mergers. 

Many studies of colour gradients try to clarify the dependence 
on the star formation history of galaxies. This aspect is of particular 
interest since it could help u s to unders t and h ow galaxies form and 
evolve. The recent study by ISuh et all j201Ch of early-type galax- 
ies drawn from the SDSS DR6 reports a tight correlation between 
the existence of steep colour gradients and ongoing residual star 
formation. This study of galaxies at z < 0.06 suggests a relation 
between elliptical gal axies with bluer core s presenting global bluer 
colours than average. iFerreras et al. found the same relation 

for spheroidal g alaxies observed by the HST at 0.4 < z < 1.5. 
iLee et al] |2008) found that steeper colour gradients appear within 
star forming galaxies, in both lat e and early-types. A similar result 
was found in the recent study by Tort ora et al.l d2010l) . where older 
galaxies are observed to have systematically shallower age and 
metallicity gradients, and thus, colour gradients. The colour gra- 



dients of spiral galaxies appear to be dominated by the fact that, on 
average, their bulges are redder than their disks. Howe ver, this is an 
overs implification for many late-type galaxies (see e.g. lBakos et al.l 
120081) . 

Cosmological simulations do not yet have sufficient reso- 
lution to make predictions for colour gradients. However, semi- 
analytical models which include a bulge-disk decomposition and 
spatially resolved star formation, could provide some use ful in- 
sight for this particular problem dStringer & Bensor]|2007l) . Nev- 
ertheless, in general, gas-dynamical simulations with high enough 
resolution to follow physical processes within galaxies are needed 
in order to produce theoretical predictions for the internal char- 
acteristics of galaxies. "Monolithic collapse" models can produce 
elliptical galaxies with colour gradients that agree with observa- 
tions if a certain sc atter is allowed for the star formation efficiency 
dPipino et al]|2010l) . These monolithic models also predict steeper 
gradients for more massiv e galaxies, something that is observed 
up to M, ~ 3.5 x 10'° M H dSpolaor et al]|2009l : feawle et alj|2010t 
iTortora et alj201oT) . For more massive galaxies, the contribution of 
mergers provides a larger scatter that smears out such correlation. 
Kobavashi (2004) used chemo-dynamical simulations of merging 
galaxies to understand elliptical galaxy formation. The resultant 
galaxies had metallicity gradients consistent with observations. The 
predicted gradients depend strongly on the relative properties of the 
progenitor galaxies, i.e., on the particular formation process of a 
galaxy. 

In this paper, we study the largest data set to date, which 
allows us to make a statistically robust analysis of the internal 
colour variation within nearby galaxies. We study four luminosity- 
threshold samples, defined in section ij2] We obtain colour gradi- 
ents directly from the SDSS photometric data and we (K+e) correct 
global magnitudes using spectral energy distribution s (SEDs) gen- 
erated with pegase dFioc & Rocca-Volmerangdl 19971) . This process 
is described in section 32.11 In section jj3]we analyse the colour 
gradients distribution and we also study the possible correlations 
between the colour gradient and global galactic properties, includ- 
ing galactic types and redshift. Finally we summarise our results 
and conclusions in section ^4] 

In this study we use Petrosian magnitudes for setting the sam- 
ple limits. Apparent magnitudes are corrected from the extinction 
due to our galaxy, but not to their internal reddening. Model mag- 
nitudes are used only for global galactic colours. All magnitudes 
used in this paper are in the AB system, unless otherwise specified. 
We assume a flat cosmology, Q. K = 0, with a critical matter density 
of fi„, = 0.3 and a Hubble co nstant, H = 100 h km s~ [ Mpc~ l with 
h = 0.7 dSpergel et. al]|2003l) . 



2 THE DATA 

Our sample is drawn from the SDSS DR7 dAbazaiian etaUl2009h 
spectroscopic main galaxy targets. The SDSS DR7 has covered 
spectroscopically an area of 9380 deg 2 , providing a redshift accu- 
racy of 30 km/ s rms for the main gal axy sample. The SD SS main 
spectroscopic sample is described by Stra uss et al.l d2002l) . Galax- 
ies within this sample have r < 11.11 (Petrosian magnitude after 
correcting for galactic extinction) and a median redshift of z ~ 0.1. 

From the SDSS main spectroscopic sample, we further select 
those galaxies that, according to the SDSS pipeline flags, have a 
reliable determination of their redshift, a good photometr)Q in all 

1 By good photometry we mean here that according to the SDSS pipeline 



SDSS DR7 colour gradients 3 



g, r, i and z bands, have a good determination of their Petrosian 
radius and a complete surface brightness profile measurement. In 
order to avoid those galaxies most affected by aperture effects, a 
lower limit to the galaxy redshift has been imposed, z > 0.01. 

After applying the above restrictions we found that one of the 
galaxies appeared multiple times due to different deblending of ob- 
jects at successive scans. After taking out this galaxy our initial 
sample contains 339711 galaxies. 

As it will be detailed in j j2.3| we quantify the internal colour 
variation using the SDSS DR7 measurements of surface brightness 
averaged within circular concentric annuli. At most there are 15 
measurements, though, in most cases there are less due to the low 
surface brightness of many galaxies. The outer radii of these annuli 
are fixed and range from 0.22 to 258.38 arcsec, growing approxi- 
mately exponentialf}0. 

We make use of the stellar masses, star formation rates, some 
spectral indices and cla ssification provided in the SPSS MPA/ JHL0 
collab oration database (Kauffmann et aL 2003a; Brinchmann et al. 
120041) , which has been updated for the SPSS PR7. These param- 
eters are obtained by extrapolating to the full galaxy the spectro- 
scopic information collected within the 3 arcsec diameter fibre set 
around the galactic centres. The extrapolation is done taking into 
account the global galactic colour and the effects due to dust. This 
extrapolation does not consider the colour gradients within galax- 
ies, which could include a systematic effect on the determination of 
global properties. 

2.1 (K+e)-correction 

The (K-t-e)-correction accounts for two effects. One, the K- 
correction, is the shift of the filter in the g alaxy's rest frame wit h 
redshift due to observing with a fixed band ( Oke & Sandage fl968h . 
The other, the e-correction, accounts for the change in luminos- 
ity due to ageing of existing stellar populations or new star forma- 
tion that galaxies experience up to a reference redshift, generally 
z = as in our case, from the epoch at which they are observed 
(Poggianti 1997). Both effects vary between the different bands. 
Within the S PSS bands we obtain larger corrections for bluer bands 
(Fuku gita et al.ll995l) . In order to make sure that we are comparing 
similar spectral regions from different galaxies we need to apply 
the (K+e)-correction. 

The (K+e)-correction is evaluated fitting the galaxy photome- 
try to the best model within a gri d of SEPs generated with the pop - 
ulation synthesis models pegase dFioc & Rocca-VolmerangeM 1997). 
The grid is made using the default template SEPs produced within 
pegase, matching the colours observed for different Hubble galac- 
tic types at z ~ 0, and includes a galaxy formed in a 1 Gyr long 
burst, an elliptical and a lenticular galaxies, spiral galaxies: Sa, Sb, 
Sbc, Sc, Sd and three irregular galaxies with different formation 
redshifts. 

There are six galaxies in our initial sample with colours that 
are not well fitted by the synthetic ones. These are left out of the 

flags the galaxies are detected at more than 5<r level, with all their pixels 
checked for possible flux peaks, not being a blended object, without too 
many interpolated pixels, without pixels hit by a cosmic ray and without a 
faulty sky subtraction. 

2 For further details we refer the reader to the SDSS help page 

http://cas.sdss. org/dr7/en/help/browser/constants.asp?n=ProfileDefs 

and the document on profiles written by R. Lupton & J. Gunn: 

http://www.astro.princeton.edu/~rhl/photomisc/profiles.ps 

J http://www.mpa-garching.mpg.de/SDSS/ . 
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Figure 1. K correction (top) and (K+e) correction (bottom) in the r-band as 
a function of redshift estimated by fitting pegase stellar population models 
to SDSS DR7 galaxy photometry. Contour levels represent an increase of a 
factor often and they are plot over single galaxy data points. 

sample. Since these objects represent 0.002% of the sample we do 
not apply any correction for their omission. 

Both K and (K+e) corrections are shown in Fig.[T]for all galax- 
ies in our initial sample. These corrections are applied to the global 
colours of galaxies in our sample. The K-correction alone tends to 
be larger for higher redshifts, but the dispersion increases at the 
same time. This behaviour has been fou nd in other works (see e.g. 
iBlanton et al.ll2003l ; lFukugita et al J 19951) . 

We can observe in Fig. [T] that the (K+e)-correction shows a 
weaker dependency on redshift th an the K-correctio n. This was also 
found in previous studies, such as lPoggiantH dl997h . The reason for 
this flattening is that as galaxies evolve, their luminosity decreases 
until a new burst of star formation occurs, being superimposed to 
an already old stellar population. 

We checked the possibility of calculating the (K+e)-correction 
individually for different galactic radii. For this purpose, we ob- 
tained magnitudes in circular concentric annuli for each galaxy, 
from the averaged surface brightness measured by the SPSS 
pipelines. However, colours in individual annuli are noisier than 
global ones and the fitting procedure basically fits noisy features 
that result in obtaining unrealistic SEPs for each annulus. There- 
fore, we obtain the global (K+e)-correction for each galaxy, which 
presumably is an average of the correction to different annuli. The 
use of an average (K+e)-correction could introduce spurious colour 
gradients. These are difficult to evaluate, since they depend on 
the intrinsic metailicity and age gradients of each galaxy, higher 
signal-to-noise data would be needed to properly compute these 
corrections as a function of radius and evaluate the effect of ob- 
serving different radii at a given projected annulus. Applying a 
global corrections implies that the resultant colour gradients would 
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Name Absolute magnitude Redshift No. galaxies 



SI 9 


M r 


- 5logh < 


-19. 


0.01 < z < 0.072 


48728 


S20 


M r 


- 5logh < 


-20. 


0.01 <z< 0.111 


82838 


S20.5 


M r 


- 5logh < 


-20.5 


0.01 < z < 0.1375 


86818 


S21 


Mr 


- 51ogh < 


-21. 


0.01 <z< 0.169 


66177 



Table 1. The given name, absolute magnitude range, redshift range and final 
number of galaxies to be studied within the defined luminosity-threshold 
samples. 



not be affected, since we are only moving the zero point on the 
colour vs. radius relation. This implies that similar colour gradi- 
ents will be found either applying the (K-t-e)-correction or only 
the K-correction. However, applying one or the other changes very 
slightly the defined luminosity-threshold samples, and has a defi- 
nite impact on the global colours associated to each galaxy. 

We will restrict our study to galaxies seen up to such a redshift 
that the (K-t-e)-corrections would not be large. Rather than the ac- 
tual value of the (K+e)-correction, the problem resides in the higher 
spread among template SEDs found at higher redshifts, which can 
introduce larger errors in the estimated correction. For the major- 
ity of galaxies up to z < 0.17, the (K+e)-corrections is well below 
a value of 0.6 in the g, r, i and z-bands. For the same range, the 
K-correction stays below 0.8 in all bands. There is another reason 
for not considering galaxies at the highest redshifts covered by the 
survey. While the distribution of the K-correction is smooth, it is 
clear that our estimate of the (K+e)-correction start to be almost 
discretised for galaxies with z > 0.25. This effect is not present at 
Z < 0.17. In order to correct galaxies at z > 0. 17 we need to enhance 
our grid of model SEDs, including different formation redshifts for 
not only Irregular galaxies but for all the Hubble types included 
as templates. This change should be done carefully, in order to re- 
produce the observed fluxes of galaxies at different evolutionary 
stages. 

Taking into account the two mentioned sources of uncertainty 
we will restrict our study to galaxies with z < 0.17. 

2.2 The luminosity-threshold samples 



We have divided our initial sample of galaxies with z < 0.17 
into four luminosity-threshold samples. These samples, though 
well defined, are not volume limited since they miss a cut in red- 
shift corresponding to the bright limit of the survey. Fig. [2] shows 
the boundaries of each sample superimposed on the distribution 
of the absolute magnitude versus redshift for all galaxies. Absolute 
magnitudes have been obtained from the (K+e)-corrected Petrosian 
magnitudes in the r-band. Table[T]summarises the characteristics of 
the luminosity-threshold samples. We will refer to each sample by 
their absolute magnitude cut, as indicated in the first column of Ta- 
bleffl 

We find that some of the galaxies within the S21 sample have 
a half light radius R50 < 0.7 arcsec, i.e., inside the first two annuli 
with surface brightness measurements. As explained in 92,3l these 
two annuli are not considered in the calculation of the colour gra- 
dients. Thus, for those galaxies within the S21 sample we will be 
calculating colour gradients with only the outer half of the galactic 
light. For this reason, we will concentrate on the other three sam- 
ples at z < 0.14, and results from the S21 sample will be considered 
with caution. 
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z 

Figure 2. Redshift versus absolute magnitude in r band for the selected 
SDSS DR7 galaxies. The different luminosity-threshold samples are en- 
closed by solid lines for the S19 sample, dashed for the S20, dotted for 
the S20.5 and dash-dotted for the S21. 

As explained in the next section 32.31 to obtain colour gra- 
dients, galaxies should have their surface brightness measured in 
at least four annuli, without considering the two innermost ones. 
Among our samples we find two galaxies which have their surface 
brightness measured in too few annuli, one within the sample 520.5 
and the other within the S19 one. These two galaxies lack a colour 
gradient measurement and are not considered further. 

As mentioned previously, we make use of some of the quan- 
tities derived from the SDSS database by the MPA/JHU collabo- 
ration. The cross identification using the position coordinates of 
the galaxies is not perfect. Fewer than 5% of the galaxies in our 
initial samples do not appear to have properties measured by the 
MPA/JHU collaboration. Only when using the MPA/JHU derived 
parameters we will be neglecting those galaxies. 



2.3 Quantifying the internal colour variation through colour 
gradients 

Different ways to compute the internal colour variation can 
be found in the literature. Some examples are: the slope of 
colour vs. normalised radius or its log arithm (e.g. Hinkl ev fe Iml 
1200 ll ; lMichardl2005l iTortora et alj|20ld) . the difference between 
colours outside and insi de certain radius (e.g. iMenanteau et al.l 
I2005L iPark & Choill2005h . the comparison between model and fi- 
bre colours when working with galaxies from the SDSS (e.g. 
iRoche et al.ll2009t) a nd the colour variation averaged to all the pix- 
els in an image (e.g. lMenanteau et alj|200lh . 

In the early stages of this work we used the SDSS DR2 to es- 
timate the best way to quantify the internal colour variation. We 
studied the internal colour variation within a sample of galaxies 
with a concentration index C > 2.63 and z < 0. 1 . The concentration 
index is defined as the ratio between two radii enclosing different 
amounts of light flux. The SDSS definition uses the radii con t aining 
90 %, R90, and 50%, R 50, of the total galactic light. lStratev3fe00lb 
and lShimasaku] l l200lh found that C = 2.6 is the boundary that sep- 
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arates early from late-type galaxies, according to both their mor- 
phological and spectral features. The actual value of this boundary 
can vary depending on level of restriction needed for contaminants 
within one or the other type of galaxies 

After thoroughly testing various methods using the quantities 
provided by the SDSS database, we reached the conclusion that the 
surface brightness profile appears to be the best way to quantify the 
internal colour variation in galaxies as opposed to magnitude based 
estimations, e.g. comparing the fibre or PSF colours with the global 
galactic colours. 

The SDSS provides measurements of surface brightness av- 
eraged within circular concentric annuli at different galactic radii. 
We obtain colour gradients from these measurements. The two in- 
nermost annuli, R < 0.7 arcsec, are excluded from the calculation, 
since they are the most affected by seeing (the median PSF width in 
the r band f or the SDSS DR7 is 1.4 arcsec). Previous studies, such 
as that from lTamura & Ohtal J2003h have reported that the internal 
colour variation is diminished if the annuli most affected by seeing 
are considered for obtaining the surface brightness profile. 

From the SDSS averaged surface brightness profile we get 
the cumulative one. This last one is fitted to a taut splinfl using 
a fourth order polynomial. Thus we need galaxies to have mea- 
surements of average surface brightness in at least four annuli. The 
interpolated surface brightness points are differentiated. We then 
convert the resulting surface brightness into magnitudes. We do this 
last step only for the radii corresponding to the outer ones from 
the initial annuli. These magnitudes are corrected with the global 
(K+e)-correction obtained for the whole galaxy. Colour gradients 
are finally obtained by calculating the slope of the straight line that 
best fits the colours, - m 2 ), at different galactic radii (those for 
which a surface brightness profile is provided in the SDSS) versus 
the radii of the annuli normalised by the galaxy half light radius in 
the r-band, R50: 



A(mi - m*i) 
MR/R50) 



(1) 



The normalisation in radius is done in order to allow us 
to compare galaxies with different light distributions and sizes. 
The colour gradients of all the galaxies within the four studied 
luminosity-threshold samples have been obtained from the sur- 
face brightness profile out to the annulus containing the radius 
R = 2 x R50. Thus, the radial coverage is consistent with previ- 
ous studies. A further analysis on the effect of radial coverage is 
made in lA2l In order to allow a fair comparison with recent liter- 
ature we have also obtained colour gradients from the variation of 
internal colour with the logarithm of the normalised radius (colour 
per radial decade). 

When computing the errors on the measured gradients we use 
a Monte Carlo method to take into account the observational er- 
rors in both the surface brightness and R50 measurements. Surface 
brightness values within a galaxy are related, with one value in- 
fluencing the values at outer radii. In order to simplify the calcu- 
lation we use a compromise by effectively considering the surface 
brightness values as independent and restricting the random values, 
obtained with the Monte Carlo method, within lcr. Thus, in order 
to include the observational errors, we repeat the whole process, 
except for the calculation of the (K+e)-correction, 1000 times for 
each galaxy using each time a random value extracted from the er- 



4 A taut spline is a piecewise cubic interpolan t for which the first and last 
internal knots are not used. See de Bood il978l) for a full description. 



ror distribution for the variables within lcr of the measurements, 
assuming that the initial errors have a Gaussian distribution. 

We obtain the probability distribution from the Monte Carlo 
runs making use of a simple estimate)^ and then we calculate the 
expected value of the colour gradient and its variance. 



3 RESULTS 

Here we explore the origin of the colour gradient distribution in 
connection with the different classifications of galaxies into either 
early and late types or star forming, passively evolving and hosting 
nuclear activity. 

We have calculated the gradients for the (g - i), (g - r), (r - i), 
(i - z) and (r - z) colours. Fig.[3]shows the distribution of the (g - i) 
and (r - z) colour gradients for galaxies in sample S20.5. We have 
found that the distributions of the colour gradients in the (g - f) and 
(g - i) colours are similar. The distributions of the colour gradients 
in the (r - i), (r - z) and (i - z) colours are also similar between 
them. As described in the introduction this is due to the similarity 
in the spectral characteristics that these colours probe. 

Fig. [4] presents the dependency of the variance in the colour 
gradient with the actual values of the (g - (') and (r - z) colour gra- 
dients. From this figure we learn that galaxies with redder cores 
(negative colour gradient) measured from the (r - z) colour gradi- 
ents, tend to have larger errors, due to the z band being the noisiest 
among the four in this study. The same is not true for the (g - i) 
colour gradients. It can be seen in Fig.|4]that steeper gradients tend 
to have larger errors, for both colour gradients. 

Fig. [3] shows that both the (g - i) and (r - z) colour gradi- 
ents for sample S20.5 have median values that are similar: -0.048, 
-0.067, respectively (-0.187, -0.277, when defining colour gradi- 
ents as colour per decade in normalised radius, in decimal logarith- 
mic units). However, the distribution for the (r - z) colour gradient 
presents a clear longer negative tail. This slight difference is likely 
to be due to the fact that redder stars contribute more in the (r - z) 
than in the (g - i) colour, giving rise to steeper gradients. However, 
the tilted distribution of the (r—z) colour gradient variance shown in 
Fig-El would also be responsible, at least in part, for this difference. 

We have also compared the distribution of colour gradients 
for the different luminosity-threshold selections defined in Table[T] 
There is not a clear change in the shape of the distributions for the 
different samples, just different numbers of galaxies. The change in 
median values from one sample to another is rather small: -0.050, 
-0.049, -0.044 for the median (g - i) colour gradient of samples 
SI 9, S20 and S21, respectively. When defining colour gradients 
as colour per decade in normalised radius we obtain the follow- 
ing median values -0.185, -0.194, -0.167 for samples S19, S20 
and S21 , respectively. These last values are roughly in agreement 
with the lTortora et alj d2010h study, where they found a non mono- 
tonic dependency of colour gradients with magnitude, when galax- 
ies fainter than M r = -20 were considered. 

In appendix [A] we analyse some of the systematic effects that 



5 The simple estimator, f, is basically a weighted histogram without a de- 
pendency on the selection of the first bin, though it maintains a dependency 
on the step size, h: 



Here the weight, w, is 0.5 if \x - x-,\ < 1 and otherwise. 
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Figure 3. Distribution of (g - i) (solid line) and (r - z) (dot-dashed line) 
colour gradients for galaxies within the sample S20.5. Median colour gradi- 
ent values are shown as short lines of the same type as their corresponding 
distribution. Histograms are normalised to give unit area underneath. 

could affect the distribution of colour gradients, including the ellip- 
ticity, the amount of light considered in the calculation and the ef- 
fect of having close pairs of galaxies. Considering the result found 
in appendix |A] and the similar distributions for all the luminosity- 
threshold samples considered in this work, we can conclude that 
the results for the colour gradient distribution are robust. 

Returning to Fig. [3] we can clearly see there that galaxies 
tend to have small internal colour variation, resulting in quite flat 
colour gradients. Nevertheless, the median of the distribution is 
negative, implying that, in general, galaxy centres are redder than 
their outer parts. This agrees with previous studi e s using smaller 
samples of galaxies (see e.g. JPeletier et al.1 Il990l : iMicharj 120051 ; 
lLeeetalf 2008; iTortora et al.ll20ld) . As described in the introduc- 
tion, Fig. [3] shows then that galaxies tend to have more metal 
rich cores that might also be hosting older stellar populations. 
Thoug h, age gradients present larger dispersion than metallicity 
ones jMacArthur et alj|20"09l ; IVlaiic et al.|[2009l ; iRawle et alj|20ld : 
lTortoraetal1l2010h . 

A remarkable feature of Fig.[3]is that the colour gradient distri- 
butions presents a single peak. This result could be indicating that 
galaxies with different structural characteristics, late and early type 
galaxies, have undergone similar star formation histories. However, 
the distribution of colour gradients presented in Fig. [3] is much 
broader than a Gaussian. We have checked whether this is due to 
the observational errors using Monte Carlo simulations. We have 
found that taking into account the colour gradients variance, we 
still obtain a distribution narrower than the observed one. Thus, we 
demonstrate that the width and shape of the colour gradient distri- 
bution cannot be explained by the error distribution alone and that 
galaxies do indeed present different intrinsic colour gradients. We 
find the colour gradient distribution to be best described as the su- 
perposition of the distributions of two populations of galaxies. 

We will explore further this last aspect in 33.21 and 93.31 Be- 
fore, we analyse the potentially different nature of those galaxies 
with extreme positive and negative colour gradients. 




Figure 4. Distribution of colour gradients errors for (g - i), (left) and (/ — z) 
(right) for galaxies in the S20.5 sample. Each contour represents a factor of 
two change in density. 



3.1 Positive and negative colour gradients 

Here we study whether extremely steep colour gradients occur only 
in certain types of galaxies. 

In Fig. [5] and Fig. [6] we explore the possibility for galaxies 
with positive and negative colour gradients to separate in a colour- 
magnitude or colour-concentration parameter space. In both plots, 
density contours correspond to the behaviour of three subsamples 
separated by their (g - i) colour gradient values. We have set as 
boundaries the 2 a values corresponding to the (g - i) colour gra- 
dient distributions of each sample. In this way we separate those 
galaxies with extreme colour gradients from the global trend. 

Fig. H] shows the absolute magnitude versus the global (g - r) 
colour, while Fig. [6] presents (u - r) v ersus light concent ration. 
Both parame t er spa ces were explored by Blant on et alj d2003h and 
iDriver et alj ( I2006I) . respectively. These studies showed the bi- 
modality of gala xies, which separate in these spaces into early and 
late-type. In the lBlanton et alj d2003h observations, the bimodality 
of SDSS galaxies is clear for galaxies fainter than M 0A 20 (ab- 
solute magnitude at z = 0.1). A similar trend is found here for the 
faintest sample 579. The brighter samples, 520, S20.5, S21, do not 
show a clear bimodality due to the magnitude cut. 

In Fig. [5] those galaxies with very steep colour gradients only 
appear at the faint end, corresponding to the region with the largest 
dispersion, but also with the largest errors in magnitude. Galaxies 
with very steep positive and negative colour are distributed simi- 
larly in Fig. [5] 

Fig. [6] shows that galaxies with extreme colours neither sep- 
arate in this parameter space. The same is seen when defining the 
colour gradient as colour per decade in normalised radius. Nev- 
ertheless, Fig. [6] shows a remarkable behaviour: the galaxies with 
the steepest colour gradients are not galaxies with high concentra- 
tion indices. Galaxies with high concentration indices are generally 
passively evolving with minimal content of gas. Since few galax- 
ies with high concentration index show steep colour gradients, then 
one expects that these steep colour gradients are more likely related 
to star formation. We will further explore this in 93.21 

From both Fig.[5]and Fig.|6]it is clear that galaxies with steep 
colour gradients, independently of their sign, appear in larger num- 
bers towards the blue end of both the (u - r) and (g - r) global 
colour distributions. Both figures show that redder and bluer cores 
appear within galaxies with similar global colours. However, for 
brighter samples, we find clearly that redder cores appear, on av- 
erage, within redder galaxies than those with bluer cores. In par- 
ticular, for the sample 520.5 those galaxies with extremely blue 
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Figure 5. The upper panel shows the absolute magnitude in the r-band ver- 
sus the global colour (g - r) for galaxies within the sample S19. Contours 
increase by factors of two in density, illustrating the behaviour of galaxies 
within different ranges of the (g — i) colour gradient distribution. Dotted 
contours correspond to galaxies with (g - i) colour gradients within the 2 <r 
limits of the total distribution. Solid/Dashed contours present the distribu- 
tion of those galaxies with their (g - ;') colour gradients above/below the 2 cr 
limit of the distribution, i.e., with bluer/redder cores than the average trend. 
The lower panel shows, in the corresponding line types, the medians and the 
distributions of the global (g - r) colour for the three subsamples selected 
according to the (g - i) colour gradients values. Histograms are normalised 
to give unit area underneath. 



(red) cores have median global colours: (g-r) = 0.82 (0.85), 
(u-r) = 2.3 (2.4), while those galaxies with flatter colour gra- 
dients present the following median global colours: (g - r) = 0.94, 
(u -r) = 2.8. 

We have also explored whether galaxies with very steep gra- 
dients appear separated when looking to parameter spaces that de- 
pend on age or metallicity against concentration index. In particular 
we have studied the variation with light concentration of the D4000 
spectral index and the Lick index Hp, both of which are mostly 
sensitive to age. We have also compared concentration to the Lick 
indices Mg2 and Fe5335, both closely related to the global metal- 
licity of the galaxy. 

We observe that both Fe5335 and Hp present values that do 
not appreciably change with light concentration, which prevents 
the separation of galaxies in such parameter space. We find a simi- 
lar bmiodality_jn_the parame ter space £>4000 vs. C, to that reported 
bv lKauffmann et all {2003b) for the SDSS DR2. Moreover, we find 
that galaxies also separate into early and late-types in the parameter 
space Mg2 vs. C. We find again that galaxies with very steep gradi- 
ents appear to have a wide range of indices values, basically cover- 
ing the same range as the global trend. It is also interesting to note 
that this is independent of the steep colour gradients being positive 
or negative. The lack of a clear locus occupied by those galaxies 
with steep gradients in the different studied parameter spaces sug- 
gest that the existence of extreme colour gradients depend weakly 
on the global age or metallicity of galaxies. 

The studied indices are measured within the SDSS DR7 fi- 
bre spectra, with a diameter of 3 arcsec. We have explored how 
our results change depending on the amount of the galactic light 



Figure 6. Similar to Fig.|5]but for (u - r) vs concentration index. 

that those 3 arcsec comprise. For this purpose, we have studied 
the behaviour of a subsample of galaxies with R5Q < 1.7 arcsec, 
i.e. galaxies with more than half of their light contained within 
the SDSS fibre spectra. This subsample exhibits similar tendencies 
to the global ones, except that in this case there are galaxies with 
steep colour gradients that have their light highly concentrated. A 
subsample of galaxies with only their cores measured within the 
SDSS spectra, with R5Q > 4.6 arcsec, give a better match to the 
global trends, in the sense that the light concentration distributes 
in a similar way to the total sample. We also find similar trends 
for the ZJ4000 vs. C, when using the spectral index derived by the 
MPA/JHU collaboration. Therefore, the tendencies of the colour 
gradients with the SDSS age and metallicity indices can be consid- 
ered as representative of the tendencies with the global ages and 
metallicities of galaxies without making any further consideration. 

We find that extreme internal colour variations tend to happen 
within galaxies that are blue, with low concentration indices and in 
the faint end of the sample distribution. Though, the difference be- 
tween galaxies with extreme positive and negative colour gradients 
is small. This result implies that if the internal colour variation was 
mainly due to supernova winds or localised star formation regions, 
then the spatial probability for these phenomena to happen should 
be quite flat within g alaxies. In their work on early-type galax- 
ies |RawleeraI] ((2010) suggested that for galaxies with a velocity 
dispersion below 140 km/ s, metallicity gradients are tightly corre- 
lated with stellar mass, while above this threshold mergers are the 
main driver for setting the colour distribution within a galaxy. Thus, 
rather than a similar origin for blue and red cores, we could be 
looking to different phenomena: AGN or gas poor mergers vs. gas 
rich mergers, occurring with a similar probability. All this points 
out the difficulty on disentangling how very steep gradients were 
originated. Further work will be needed to understand the similar 
distributions of positive and negative colour gradients. 

3.2 Colour gradients within early and late type galaxies 

In this section we study the internal colour variation within early 
and late-type galaxies, using different parameters to classify galax- 
ies into one of those categories. 
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Figure 7. Variation with the concentration index (left), (g - r) colour (middle) and specific star formation (right) of the (g - i) colour gradient, top panels, and 
its error, bottom panels, for galaxies within the S20.5 sample. The panels showing colour gradients have superimposed median colour gradients values and 
their lcr dispersion range. The panels showing the variance of the colour gradients present density contours with increments of a factor of two. 



The left panel in Fig.UJshows the variation of the (g - ;') colour 
gradient with the concentration index, C, for galaxies in the 520.5 
sample. As it can be seen in this plot, median colour gradients tend 
to be flatter for galaxies with higher values of the light concentra- 
tion. In fact, the same is seen for all the studied colour gradients: 
(g - i), (g - r), (r - z), (r - i), (i - z), in the different luminosity- 
threshold samples considered here. This was expected given the 
trends seen in the previous section. Due to the dispersion, this cor- 
relation is not statistically robust, with a Spearman coefficient of 
0.19. However, the increase in the dispersion of observed colour 
gradients is not correlated with a trend in the gradient errors and 
so, it is significant. 

We find that early-type galaxies tend to have flatter colour gra- 
dients than late-type galaxies. Most spirals are classified as late- 
type and have an exponential light profile. In simple terms, these 
galaxies are formed by a bulge and a disk component. Bulges share 
features with early-type galaxies, being dominated by an old stel- 
lar population. Meanwhile, disks contain, in general, regions of star 
formation, providing an increase in the percentage of young stellar 
populations. These differences allow for greater colour variations 
to happen in comparison with early-type galaxies, whose charac- 
terist ics are much more homogeneous. 

iTortora et alj 1 120101) also found a similar trend for colour 
gradients with the Sersic index, which is tightly related to the 
concentration index. T hese results also agree with the study by 
IWeinmann et alj d2009h . were they found that galaxies with C < 3 
present steeper colour gradients. 

Colour has also been widely used in the literature to split 
galaxies into early and late-types. The middle panel in Fig.[7]shows 
the variation of the (g - i) colour gradient with the global (g - r) 
colour of galaxies. In agreement with Fig.[5]and Fig.[6] it is clear in 
Fig.[7Jthat the redder a galaxy is the flatter the colour gradient is ex- 
pected to be. The trend is similar to that seen for the concentration 
index: early-type galaxies present flatter colour gradients. However, 
the correlation with global colours is weaker than with the concen- 
tration index, with Spearman index of ~ 0.10. This could indicated 
that colour gradients depend more strongly on morphological char- 
acteristics than on global colours. For example, if the colour gradi- 
ent is set by the relative contribution of bulge and disc for late-type 
galaxies, one would expect that it would correlate more strongly 
with concentration index than with global colour. 

For a given mass range, late-type galaxies are, in general, 
forming stars more efficiently than early-types ones. The right 
panel in Fig. [7j presents the variation of colour gradient with the 
specific star formation derived by the MPA/JHU collaboration. 



With a Spearman index of ~ 0.14, the correlation between these 
two parameters is not statistically robust, likely due to the increase 
of the dispersion for those galaxies with higher specific star forma- 
tion rates. The variance of the colour gradients does not determine 
the dispersion increase, which makes, again, the tendency signifi- 
cant. 

Assuming a constant star formation, the inverse of the specific 
star formation, M, IS FR, will provide an estimation of the time it 
took a galaxy to build up its observed stellar mass. Thus, higher 
values of specific star formation will correspond to younger galax- 
ies, or at least, with a larger percentage of young stars. The found 
correlation between the specific star formation and the colour gradi- 
ent indicates that steeper gradients appear for galaxies with higher 
percentages of young stars. 

From Fig. [7] we can conclude that early-type galaxies have 
smaller internal colour variations compared to late-type galaxies. 

We have also compared the shape of the distribution of colour 
gradients within early and late-type galaxies, selected accordingly 
to either concentration index, global colour or specific star forma- 
tion, finding that early-type colour gradients have narrower distri- 
butions than late-types. As can be seen in Fig. [8] this is also true 
when separating galaxies into early and late-types using other pa- 
rameters. 

Fig. [8] shows the (g - f) colour gradient distribution for galax- 
ies from t he S20.5 sample separated into early and late-types fol- 
lowing the Kauffmann et alj j2003bh work, by using the following 
parameters derived by the MPA/JHU collaboration: the £)„4000 
index, sampling the spectral 4000A break, the Balmer index H 6A 
and the logarithmic stellar mass. Mass does not clearly sepa- 
rate galaxies into early and late-types, but the characteristics of 
galaxies with masses bel ow or above 10 105 M o are rather different 
i lKauffmann et alj[2 003b). Despite median values being quite close 
for early and late-types galaxies, it is apparent that late-type galax- 
ies present a wider distribution with more extended tails than for 
the case of the early-typ e galaxies distributio n. 

As it is found here. lTortora et alj (2010) also report ed a larger 
disper sion for colour gradients among late-type galaxies. ISuh et all 
d2010h found that although most of the early-type galaxies have 
quite flat colour gradients, there is ~ 30% with residual star forma- 
tion that produces steeper colour gradients, in consonance with the 
distributions reported here. 

Following 93.11 we have also studied here the relation between 
global galactic colours and the existence of extreme colour gradi- 
ents within galaxies with low and high concentration indices. We 
find again that steep colour gradients appear within bluer late-type 
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Figure 8. Distribution of the (g - ;') colour gradients within galaxies from 
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underneath. 



galaxies. However, in the case of early-type galaxies the global 
galactic colour does not appear to be correlated with the existence 
of extreme colour gradients. 

When defining colour gradients as the variation of colour with 
the logarithm of the normalised radius, we obtain that early-type 
galaxies (C > 2.6) have a median (g - i) colour gradient of -0.13 
for both the S19 and S20.5 samples, while late-type galaxies have 
median values of -0.27 and -0.32, for those two samples re- 
spectively. The quote d values for earl y -type galaxies are smaller 
than those found by La Barbera et al] J2010h . In the study from 
La Bar bera et aDHoiO), they apply two cuts in order to minimise 
the contamination of the early-type galaxy sample. Since we have 
only separated galaxies on basis of their concentration index values, 
we could have a larger fraction of late-type conta minants, providing 
steepe r colour gradients than those reported by lLa Barbera et al.l 

In summary, the probability for finding strong internal colour 
variations within early-type galaxies is lower than for late-type 
galaxies. 



3.2.1 Early and late type galaxies based on their colour 
gradients 

IChoi et alj J2007l) have studied the behaviour of colour gradients 
with absolute magnitude for early and late-type galaxies, separated 
by their position i n the colour gradie nt versus colour parameter 
space as defined bvlPark & Choil (2005). In order to compare our re- 
sults with the lChoi et al.l j2007l) study we have also separated galax- 
ies into early and late-type according to the IPark & Choil (2005) 
method. Our definition of colour gradients is different since they 
just compare the light outside and inside half the Petrosian radius 
of galaxies. Moreover, they use magnitudes K-corrected to redshift 
0.1. We have not applied any corrections to take into account these 
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Figure 9. Both panels show data from galaxies within the S19 sample. 
Galaxies have been separated into e arly, in black, and late-type galaxies, in 
grey, following IPark & C hoi (2005) criterion. Contours represent a factor 
of two change in density. Left: (u — r) colour vs. (g - i) colour gradient. 
Blue contours show the density for galaxies with C < 2.6, those in red 
correspond to galaxies with C > 3. Right: Absolute magnitude in r band vs. 
(g - r) colour gradients. 



differences, we will thus compare tendencies qualitatively, rather 
than actual values. 

The left panel in Fig. [9] shows the (u - r) colour versus the 
(g - i) colour gradient, using black/grey dots for those early/late- 
type galaxies according to the lPark & Choil d2005h criterion. In the 
same panel we have superimposed density contours for galaxies 
with different lig ht concentrat i ons: C < 2.6 and C > 3. We can 
see there that the Par k & Choil {2005) classification scheme agrees 
quite well with separations based on the concentration index. 

The right panel in Fig.|9]shows the (g - i) colour gradients as 
a function of absolute magnitude. Fig. [9] shows that those galax- 
ies classified as early-types have a very constant median colour 
gra dient with absolute mag nitude. The same behav iour is reported 
by IChoi et alj to07l) and lLa Barbera et al.1 to05t) thoug h other 
studies find some trends (see e.g. lLa Barbera et al. I2010Q . Early- 
type galaxies within the 579 sample, especially those fainter than 
M r 21, present values with increasing dispersion towards pos- 
itive values, i.e. towards having bluer cores. This could be related 
to the branch of early-type galax ies with residual sta r formation, 
as the E+A galaxiefl selected by IPark & Choil {2005) method. In- 
deed, in their work they showed that one of the advantages of their 
way of separating galaxies is that E+A galaxies are classified as 
early-types, although they are bluer and less concentrated than most 
early-type galaxies. In fact, this increase in dispersion towards bluer 
cores within early-types is the only difference found in the colour 
gradient vs. absolute magnitude parameter space when using the 
concentration index to separate galaxies into early and late-types. 

Unlike global galactic colours, the internal colour variation 
do not relate clearly to the luminosity of galaxies, though median 
colour gradients tend to be steeper at brighter absolute magnitude 
bins. For fainter magnitudes than the ones studied here. lChoi et al.1 



6 E+A galaxies are dominated by an old stellar population but their spec- 
tra show signatures of the occurrence of minor s tar formation bursts that 
have recently ceased. See lYam auchi & Goto 1 20051) for a study of colour 
gradients in E+A galaxies. 
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Figure 10. Variation with stellar mass of the median (g — i) colour gradient 
for early, circles, and late-type galaxies, stars, within sample S19. Error 
bars corresponds to the lcr range in dispersion. The bottom panel shows 
the variance of colour gradients within early, solid contours, and late-type 
galaxies, dotted contours. Plotted contours increase by a factor of two in 
density. 



j2007h observed that the colour gradients of late-type galaxies in- 
crease with magnitude. Nevertheless, within the magni tude range 
studie d here the median colour gradient s reported by Choi et ail 
(2007) also stay quite flat with magnitude. lBakos et alTd2008h stacl 
ied the internal colour variation within 85 spiral galaxies finding 
that it is almost constant for a given luminosity. 

In comparison to early-types, late-type galaxies have steeper 
colour gradients and present a larger dispersion that in creases for 
fainter magnitudes. These results are in agreement with lChoi et al.l 
d2007h and in consonance with our previous findings. 



3.2.2 Velocity dispersion and stellar mass 

Due to the complexity for obtaining the velocity dispersion, the 
SDSS provides measurements only for galaxies with a clean spec- 
trum typical of an early-type at z < 0.4. The subsample of galax- 
ies with a reliable measurement of the velocity dispersion presents 
slightly flatter colour gradients with larger velocity dispersions. 
Such tendency is not statistically robust for the (g - i) colour gra- 
dient, with an Spearman index of ~ 0.12, due to the increase in 
dispersion. Though, as for the parameters studied before, the dis- 
tribution of errors alone cannot explain the change in dispersion, 
which makes the tendency signific ant. Th e same tendency has been 
found in the study by Roc he et all d200°h on early-type galaxies. 

Stellar mass and velocity dispersion are tightly correlated. The 
trend seen for the velocity dispersion is supported by the correla- 
tion that we find between the (g - i) colour gradient and the stellar 
mass derived by the MPA/JHU collaboration. For a similar range in 
both masses and luminosity, our result for the most massive galax- 
ies to present flatter colour gradients agrees with previous studies 
dWeinmann et al.ll2009l : [Suh et al.ll2010| - |Tortora et alj|2010h . 

Fig. [10] shows the variation with stellar mass of the median 
(g - i) colour gradient for early, with C > 2.6, and late-type galax- 
ies, with C < 2.6. The median of the colour gradient distribu- 
tions for both type of galaxies remain quite flat for masses above 



M > 10 95 M o . Less massive galaxies tend to present flatter or even 
positive colour gradients, although the dispersion is much larger. 

Our sample of galaxies contain rather massive galaxies, with 
M > 10 9 M o . For this rang e of masses there is a w ide dispersion 
for the star formation rates dBrinchmann et al]|2004l) . thus it is dif- 
ficult to connect this result with the different mechanisms that are 
able to produce colour gradients. Splitting our sample into early 
and late types we find that only those early-type galaxies with 
M, < 10 9 5 M o within the S19 samp le present a median positive 



colour gradient. Spol aor et al.l d200°l) have studied the dependency 
of metallicity gradients with stellar mass finding them to be tightly 
related if M, < 3.5 x 10 10 M a , as predicted by pure monolithic 
collapse models (e.g. IChiosi & Canard l2002h . while more mas- 
sive galaxies do no t appear to show this relation, as other studies 
previously proved (|Michard 2005; Sanchez-Blazque z et al. 1 120061: 
iTortora et alll201fj|) . For those galaxies in the lowest stellar mass 
range, the metallicity gradient becomes shallower with decreasing 
mass, which could translate into age gradients having a dominant 
effect on the colour gradients. This will provide the possibility to 
invert the colour gradients. In general, those early-type galaxies 
at z = with the lowest mass c ould present higher s tar forma- 
tion rates than more massive ones dDe Lucia et al . 2006), which in 
turn could provide the existence of steeper positive colour gradients 
(bluer cores). We find a similar result for late-type galaxies, in the 
sense that only those galaxies in the lowest stellar mass bins present 
null or positive median colour gradients values. 

The tendency fou nd for ear l y-type gala xies in Fig. 1 101 agree s 
with the studies from ISuh et al.l hoioh and ITortora et alj d2010h . 
supporting that bluer cored early-type galaxies only appear in the 
low mass bins. 

In the case of late-type galaxies, localised regions of star for- 
mation can affect the colour gradients but dust and stellar migration 
can also have a significan t impact on them dSanchez-Blazquez et al.l 
120091: IVlaiic etalj|2009l) . Stellar migration happens immediately 
after the formation of the wave density that leads to the formation 
of spiral arms, but also after the formation of bars or disk insta- 
bilities that can be triggered by tidal forces. Stellar migration is 
strongest for the most massive galaxies and tends to wipe out ini- 
tial metallicity gradients, and, therefore, colour ones. This effect 
rises age gradients opposed to the m etallicity ones, thou gh the dis- 
persion for these increase with mass dRawle et alfeOld) . The com- 
bined effect of age and metallicity gradients for massive galaxies 
can explain why we find the colour gradient to stay quit e flat for 
late-type galaxies of different masses. rTortora et aTl d2010f) found a 
clear steepening with mass for the colour gradient within late-type 
galaxies. As mention before, we find that the least massive galax- 
ies in our sample present flatter median colour gradients, i.e., we 
find a very sl i ght tre nd in the same direction as the one reported by 
ITortora et alj 1 20 101) . However, the slopes of these trends are differ- 
ent. The origin of this difference is unclear but it could be related 
with the different definitions of colour gradient or the sample selec- 
tion. 



3.3 Colour gradients, star formation and nuclear activity 

Beside stellar winds and migrations, the nuclear activity could also 
be responsible for the existenc e of anomalous internal c olour varia- 
tions within galaxies. In fact, iMenanteau et alj (2005) found one 
blue nucleated spheroid with star formation associated with the 
presence of an AGN. 

Here we try to disentangle the effects of star formation and 
nuclear activity as responsible for setting colour gradients. We 
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Figure 11. Distribution of the (g - r') colour gradients within galaxies from 
the SI 9 sample, separated into AGN plus LINER S (solid line) and star form- 
ing galaxies (dashed line) following Brin chmann et ail 120041) . The top left 
panel shows the distribution for all the galaxies in the sample, while the 
rest show the distribution within the mass range specified in the legend. 
Arrows show median values. Histograms are normalised to give unit area 
underneath. 

use the classification done by the M PA/JHU collaboration fo r 
the SDSS DR7 based on the one from iBrinchmann et at] J2004I) . 
which makes use of the location of a galaxy in the parame- 
ter space log([OIII] 5007/H B ) versus loe ([NII]6585/H a ), the so 
called BPT diagram (Bald win et al.|[l981l) . taking into account the 
signal to noise (S IN) of the spectral lines. In this parameter space 
it is possible to separate star forming galaxies from those with nu- 
clear activity, though there is a region were galaxies are considered 
as composites. Due to the requirement of a minimum S IN value for 
the lines, not all galaxies can be classified following this scheme. 

We have explored the locus occupied in the BTP diagram by 
those galaxies with very steep colour gradients, finding that they 
distribute following the mean trend. 

Fig.QT]presents the distribution of the (g - i) colour gradients 
within galaxies from the SI 9 sample, separated into galaxies with 
nuclear activity and star forming galaxies. The left upper panel of 
this figure presents the whole sample, showing that steep colour 
gradients are more likely to appear within star forming galaxies 
rather than those with nuclear activity. In fact, there is a minute in- 
crease, ~ 2%, in the percentage of galaxies classified as star form- 
ing among those with steep colour gradients, compared with the 
mean trend. This last tendency goes in the direction expected from 
the above discussions. 

The studied early-type galaxies present a higher percentage of 
galaxies with nuclear activity, independently of their colour gradi- 
ents. Thus, the top left panel in Fig.[TT]could be driven by the differ- 
ent distribution of colour gradients for early and late-type galaxies 
seen in 93.21 

The remaining panels in Fig. [TT] explore the colour 
gradient distribu t ions fo r galaxies within three mass ranges. 
Kauff mann etai] J2003cl) found that different levels of nuclear 
activity appear in galaxies with different masses. Galaxies with 
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colour gradients within galaxies in sample S20.5. 

masses below 10 l0 M o could present a weak AGN, but not a strong 
one. In this mass range we observe that those galaxies from the S19 
sample that are hosting an AGNs present a higher percentage of 
bluer cores than those forming stars. The other studied samples of 
galaxies do not present this difference, but yet, they contain fewer 
low mass galaxies. The observed trend could be explained if those 
weak AGN were triggering star formation localised only in the cen- 
tral region of their host, or, more unlikely, if the AGN within these 
low mass galaxies was driving dust out of the central galactic re- 
gion, lowering its internal extinction, providing bluer cores. 

Galaxies more massive than 10 10 M o , host most of the AGN 
activity. As expected from 93.2.21 the colour gradient distribution 
for the most massive galaxies is quite narrow, with few galax- 
ies presenting steep colour gradients. The median colour gradients 
for massive galaxies remain quite constant and comparable to the 
global median. 

Our results suggest that nuclear activity is a marginal driver 
for creating steep colour gradients, and that it is overcome by the 
dependency between colour gradients and morphological type, at 
least for those galaxies with masses above 10 10 M o . 

3.4 Colour gradient change with redshift 

Comparing the two panels in both Fig.[5]and Fig.[6]for the samples 
under study, we have observed that galaxies with steep colour gra- 
dients tend to be redder when moving to samples including galaxies 
at higher redshifts. The global sample also gets redder when going 
from sample SI 9 to S20.5, but the variation is smaller. We have 
observed that colour gradients do not depend strongly on absolute 
magnitude, except for an increase in the scatter at the faintest end. 
Thus, this result could suggest a stronger evolution for those galax- 
ies with steeper colour gradients in the recent past. 

Fig. [12] shows the variation of the (g - i) colour gradient with 
redshift for early and late-type galaxies within the 520.5 sample. 
For both types of galaxies the median colour gradients stay constant 
with the redshift, while the dispersion and errors increase. We find 
that for a fixed magnitude the dispersion in the colour gradient as 
a function of redshift is hi ghly reduced. Th e opposite is not true. 
This is in agreement with Bak os et at] d2008b . 

We observe that it is more probable to find a galaxy with a 
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steeper colour gradient at higher redshift. However, this trend of 
the dispersion could be driven by the errors distribution. 

We have explored visually 48 early-type galaxies randomly 
selected among those with positive gradients, outside the lcr dis- 
persion range, and within the median value to explore the possibil- 
ity of those early-type galaxies with steeper colour gradients to be 
more affected by artifacts. We find the same amount of galaxies ei- 
ther close to stars or other galaxies among both groups of galaxies. 
Thus, the observed tendency does not appear to be due to a higher 
presence of artifacts among those early-types with steeper colour 
gradients. 

In agreement with our study, lLa Barbera et al. I d2003h found 
that colour gradients within spheroids are weakly dependent on the 
re dshift, up to z ~ . 65. C olour gradients within the disks studied 
by lLa Barbera et al J (2003) change appreciably from z ~ to z ~ 
0.65. However, this change is not so clear within the redshift limits 
used in our study. 

The fact that median colour gradients do not appear to evolve 
at z < 0.17, points out that most of the star formation activity oc- 
curred previously. 



4 SUMMARY AND CONCLUSIONS 

In this paper we have studied the internal colour variation within 
galaxies from the SDSS DR7, with both clean photometry and 
spectroscopic measurements. We have obtained the global (K+e)- 
correction by using SED templates from pegase and applied them 
to each galaxy in order to divide our initial sample into luminosity- 
threshold ones, with redshifts up to z = 0.17. These samples are the 
largest used for studying the internal colour variation in galaxies. 

The internal colour variation has been quantified through the 
slope of the relation between radial colour and the galactic radius 
normalised to the one containing half the total galactic light, R50. 
Colour gradients have been obtained from the SDSS measurements 
of surface brightness circularly averaged in concentric annuli. Ob- 
servational errors have been taken into account using a Monte Carlo 
method in the colour gradient calculation. This approach is inde- 
pendent of previous ones in the literature, therefore the agreement 
with earlier works underscores the robustness of the colour gradi- 
ent as a useful parameter for studying particular aspects of galactic 
evolution. 

In agreement with previous studies, we have found that colour 
gradients distribute around median values that are negative but 
close to zero, i.e. galaxies, independently of their morphological 
type, tend to have slightly redder cores than their external parts. 

We have explored some of the possible systematics that could 
affect the measurement of the colour gradients. In particular we 
have explored the possible dilution due to the fact that we start our 
calculation of the colour gradient from a circularly averaged surface 
brightness profile. We observe that very elongated galaxies tend to 
have steep colour gradients, the op posite trend to the one expected 
if they were diluted. According to lChoi et al.l ( 120071) . steep colour 
gradients within elongated galaxies are due to their internal dim- 
ming. In the calculation of the colour gradients we have left out 
the innermost annuli, since they are the ones most affected by see- 
ing. We have explored the effect of obtaining colour gradients using 
different number of annuli, finding that for the range in which we 
work our results on colour gradients are robust. We find that galax- 
ies with larger half light radii tend to have steeper colour gradients, 
being an intrinsic characteristic of galaxies and not an artifact. 

We have selected a subsample of close pairs of galaxies that 



appear to be interacting. The sample is small, far less than 1% of the 
total sample, and it represents an inferior limit of such close pairs 
of galaxies with respect to what the total percentage that might be 
found when considering the whole photometric SDSS sample. De- 
spite the small number of galaxies in this sample, we find a larger 
proportion of galaxies with bluer cores here than in the whole pop- 
ulation sample. As close pairs are more likely to be interacting and 
normally these interactions are associated with an enhancement of 
star formation, this supports the connection between steep colour 
gradients and recent bursts of star formation. 

We find that the colour gradient distribution does not follow 
a Gaussian due to an excess of galaxies in the tails that cannot be 
explained alone by the errors distribution. A double-Gaussian fits 
better the colour distribution, pointing out a possible contribution 
of two distinctively different galaxy classes. Separating the initial 
sample of galaxies into early and late-types we find that the dis- 
tribution for colour gradients within late-type galaxies is more ex- 
tended than that for early-types. Nevertheless, the distribution of 
early-type galaxies is still broader than a Gaussian. We find this ex- 
cess of steep colour gradients to be likely related to age gradients, 
AGN activity or an unusual dust content. 

The observed distribution of colour gradients can be under- 
stood as a consequence of entangled effects. Gradients produced 
in mergers depend on the time scales that interplay in their cre- 
ation and the mixing of elements that can contribute to their de- 
struction. Galaxies observed at different epochs after a gravita- 
tional interaction can present colour gradients evolving with time. 
Nevertheless, even when observing galaxies at a similar time af- 
ter a merger one expects to find a distribution of colour gradients 
and not a single value, since the metallicity enrichment appear to 
be highly non-un iform, due to the multi-pha se nature of the inter 
galactic medium dSpringel & Hernouistll2003l) . Besides. iKobavashU 
( 2004) simulated the formation of early-type galaxies through dif- 
ferent kind of mergers, finding that the metallicity gradients (the 
main origin of colour gradients) strongly depend on the ratio of 
stellar masses and gas content between the progenitors. Moreover, 
even if early-type galaxies were formed in a "monolithic collapse", 
for a given stellar mass there would be certain scatter in star for- 
mati on efficiencies wh ich provides a scatter in colour gradients 
(e.g iPipino et al.l 2010). For late-type galaxies, the colour gradi- 
ent dispersion is further enhanced due to the dust content and the 
bulge to disk relation depending on the observing inclination angle. 
The complicated patterns of stars migrations that follow the spiral 
arm formation can also enhanc e, at least temporally, this dispersion 
dSanchez-Blazquez et alj|2009f) . 

Galaxies with steep gradients share some global properties in- 
dependently of presenting bluer or redder cores. Nevertheless, red- 
der cores happen, on average, within galaxies redder than those 
with bluer cores. We have also observed that galaxies with high 
concentration values do not present steep colour gradients. It is 
known that mergers increase the galactic size and that galaxies with 
highly concentrated light have features indicating either a passive 
evolution in their recent past or the lack of gas, ruling out the possi- 
ble existence of great internal colour variations. Thus, both results 
suggest the existence of a relation between colour gradients and 
the recent formation history of galaxies, in the sense that very steep 
colour gradients are likely related to recent star bursts. 

We have found that galaxies with steeper colour gradients are 
more likely to be late-types according to either their structural pa- 
rameters or to their specific star formation, i.e., a higher presence of 
young stars. It is interesting to note that we have found that colour 
gradients are less tightly related to global colours than the con- 
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centration index. IPark & Choil d2005h found out that the parameter 
space colour-colour gradient could be very useful to select galaxies 
matching their optical morphology at low redshifts. Our study indi- 
rectly supports theirs, suggesting that colour gradients origins are 
more related to the recent formation history of galaxies, the kind 
of gravitational interactions recently suffer ed, than with the globa l 
galactic colours. Different studies (see e.g. lSchawinski et"al1l2009l) 
are finding spheroidal galaxies with blue colours and red spirals. 
This proves the usefulness of a parameter like the colour gradient 
that is sampling something that cannot be seized through the global 
colour. 

We have found that the possible correlation between colour 
gradients and nuclear activity is overcome by the stronger relation 
with galactic type, since nuclear activity appears more common 
among early-types in our sample of galaxies. 

At intermediate redshift, 0.37 < z < 0.83, Menanteau et all 
d200lh found that 30% of field spheroids (classified based on their 
morphologi cal parameters) s how an unusual internal colour varia- 
tion. Menanteau et al. (2001) associated this large colour variation 
with either star formation or AGN activity, in consonance with our 
results for a closer sample of galaxies. 

Within the redshift range studied here, the colour gradients 
remain rather constant, implying that the studied galaxies had their 
last residual bursts of star formation at earlier times. 

The results presented in this work agree with previous stud- 
ies, though, here we are more dominated by the scatter due to the 
larger sample explored. Nevertheless, our results are very robust 
since they hold for the four studied luminosity-threshold samples. 

Colour gradients are normally defined in the literature as the 
colour variation as a function of the logarithm of the radius. Here 
we have used the variation of colour as a function of the radius it- 
self. Using the former definition we have checked that we obtain 
qualitatively similar trends for colour gradients, except that their 
distribution covers a broader dynamical range. For the 520.5 sam- 
ple, the median values of these distributions are -0.187 and -0.277, 
respectively for the (r - z) and (g - i) colour gradients. In our cal- 
culation of colour gradients we assume the observational errors in 
both colours and radii to be normally distributed. When using the 
logarithm of normalised radius, the error distribution deviates from 
a Gaussian. Thus, in order to make an accurate comparison we 
should modify the way we are including the observational errors 
in our calculation. This is likely to have a fundamental effect on the 
estimated errors, while colour gradients would remain practically 
unchanged. Leaving apart this issue, when defining the colour gra- 
dient as colour variation as a function of the logarithm of the radius 
we find qualitatively similar tendencies to the ones previously re- 
ported. The correlations appear stronger when defining colour gra- 
dients in this way, due to the enhanced dynamical range. 

The fact that steep colour gradients are linked with a higher 
presence of young stars could provide an alternative method for 
exploring the interactions between galaxies. A thorough spectral 
analysis of those galaxies which present extreme colour variations 
could help us to better understand the mechanisms that generate 
either bluer or redder cores within galaxies. 
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APPENDIX A: POSSIBLE SYSTEMATIC EFFECTS IN 
THE COLOUR GRADIENTS 

Al Effect of the ellipticity 

We obtain colour gradients from surface brightness measurements 
that are circularly averaged. Thus, colour gradients within elon- 
gated galaxies could be diluted. To address this issue we have sep- 
arated galaxies within sample S20.5, the most populated one, into 
different ranges according to the ratio between the projected minor 
and major axes of the galaxies in the r-band, b/a. The top panel 
of Fig. IA1 1 shows the variation of the (g - i) colour gradient with 
b/a. We can see that more elongated galaxies present steeper colour 
gradients, with median values varying from -0.07 for galaxies with 
b/a < 0.4 to -0.04 for b/a > 0.8. This tendency appears indepen- 
dently of the band used to define b/a. The same tendency of steeper 
gradient for more elongated galaxies is found for galaxies in the 
other studied samples. It is also seen when colour gradients are de- 
fined as the variation of colour per logarithmic normalised radius. 
However, this is not observed for the (r—z) colour gradients, where 

all median values are similar, 0.07, independently of their b/a 

value . 

The bottom panel of Fig. lAll shows that the distribution of the 
error does not appreciably change as a function the projected axis 
ratio. Thus, the error distribution is not enough to explain the slope 
change of the median values around b/a < 0.3 seen in the top panel 
of Fig.lATI 

This result shows that the dilution of colour gradients due to 
the use of circularly averaged surface brightness does not erase 
colour gradients although it may still dilute them somewhat. 

Galaxies in the range b/a < 0.4 are late-types. Very inclined 
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Figure Al. Variation of (g - i) colour gradients with b/a in the r-band, top 
panel, and its errors, bottom panel, for galaxies within the S20.5 sample. 
In the top panel we have superimposed median colour gradients values and 
their Icr dispersion range for b/a in the r-band, circles with solid error bars, 
but also for b/a in the g-band, triangles with dashed error bars. The bottom 
panel shows density contours with increments of a factor of two. 

late-type galaxies are most affected by their internal absorption, 
something that could explain their steeper colour gradients. Due 
to their small gas content, early-type galaxies properties are not ex- 
pect ed to be affected b y their inclination. 

Cho i et al. I d2007t) studied the dimming effect for inclined late- 
type galaxies finding that it is better to exclude those with b/a < 0.6 
when studying properties that could be affected by galactic internal 
dimming. Considering only galaxies with b/a > 0.6, we find that 
the linear variation of the median of the (g - i) colour gradient 
with b/a has a slope smaller than 0.02, for all of the luminosity- 
threshold samples considered in our study. As expected from the 
large dispersion found in Fig. IA1I the Spearman coefficient is too 
small to support statistically such a linear relation. 

We can conclude that colour gradients do not vary in a signif- 
icant manner with b/a, except for very inclined late-type galaxies, 
which are most likely affected by their internal dimming. Since the 
z-band is the least affected by dust extinction, this last result is con- 
sistent with the lack of (r - z) colour gradient variation with b/a. 



A2 Effect of the number of annuli considered 

Here we study the bias that we could be introducing by leaving 
out the innermost light, R < 0.67 arcsec, of the galaxy. For all the 
studied galaxies, colour gradients have been obtained from the sur- 
face brightness within two times the radius containing half the total 
galactic light. In section ^2]it was pointed out that the S21 sample 
contains galaxies with half light radii, R50, within the two inner- 
most annuli that we neglect computing the colour gradient. For the 
rest of the samples this does not occur. Moreover, in these sam- 
ples less than 5% of the galaxies have their half light radii within 
the first annuli considered in the calculation of colour gradients, 
0.67 < R50 (arcsec) < 1.0 (3 rd annulus in the SDSS database). 

Fig. I A2I shows the distribution of colour gradients of galaxies 
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Figure A2. Distribution of (g - i) colour gradients for galaxies with 
b/a > 0.6, within the S20 sample separated by the radius containing 
half their light: R5Q < 1.7 arcsec (outer radius of 4'' 1 annulus, solid 
line), 1.7 < i?50 (arcsec) < 3 (limits of the 5' ,! annulus, dotted line), 
3 < R50 (arcsec) < 4.6 (limits of the 6' 1 ' annulus, dashed line), R50 > 4.6 
arcsec (outside the 6' 1 ' annulus, dot-dashed line). 



with b/a > 0.6 from the 520 sample, divided in four ranges of half 
light radii. The cut in ellipticity allows us to isolate the tendency 
of colour gradients with R50. The ranges of R50 have been defined 
using as boundaries the outer radii of the 4' h , 5 ,h and 6' 1 ' annuli 
in which the SDSS database provides averaged surface brightness 
measurements. Fig. lA2l shows that the largest amount of galaxies in 
sample 520 have their R50 within the 5'' 1 surface brightness annuli. 
This is also the case when considering all galaxies, independently 
of their ellipticity and it also holds for galaxies within the sample 
S20.5, while for sample 5/9 there is an even amount of galaxies 
with their R50 in either the 5 th or 6 th annulus. As it was pointed 
out before, all galaxies in our samples have at least four surface 
brightness measurements, from the 3' rf to the 6' 1 ' annulus. Thus, 
except for the 52 / sample, when obtaining colour gradients we are 
taking into account surface brightness measurements both inside 
and outside R50. This point is essential for any robust conclusion 
about the behaviour of colour gradients. 

Median colour gradients are steeper for higher values of 
R50 (arcsec), for all the defined luminosity-thr eshold samples. The 
same tendency has been found in the study bv lRoche et alj ( 2009) 
on early-type galaxies. 

In Fig. lA2l the ranges of values we study are defined by half 
light radii values in arcsec. Therefore, gradient values could be 
smeared out depending on the amount of light considered. To fur- 
ther explore this point, we have also studied the tendency of colour 
gradients when comparing galaxies with similar physical transverse 
radius, i.e., defining ranges of R50 measured in hr l kpc. Different 
amount of surface brightness annuli account for a different amount 
of light for galaxies with similar sizes at different distances. Thus, 
in this way galaxies with colour gradients extracted from different 
amount of annuli are included in the same range. 

Independently of the galactic ellipticity, the same tendency 
found for R50(arcsec) applies when separating galaxies by their 
-R50(/r'fcpc): the larger R50 is, the steeper a colour gradient 
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Figure A3. Similar to Fig. lAll for the radius containing half the total light 
(hT l kpc) within galaxies from the S20 sample. 



is likely to be. Galaxies within the S20.5 sample have median 
(g - i) colour gradients: -0.036, -0.048, -0.061, -0.080 for 
R50(h-'kpc) < 2; 2 < R5Q(hr l kpc) < 3; 3 < R5Q(h~ l kpc) < 4; 
R50(h~' kpc) > 4, respectively. The same tendency for median val- 
ues is found for galaxies within the SI 9 and 520 samples. In Fig. [A3] 
we present the variation with R50(lr l kpc) of the (g-i) colour gradi- 
ent and its error, for galaxies in 520. Median colour gradient values 
tend to decrease with larger half light radii, though, the correlation 
is dominated by the dispersion. When defining the colour gradient 
in decades of normalised radii (radius in logarithmic units), we also 
find that colour gra dients get steeper for galaxies with larger R50. 
iTortora et alj d2010h finds a similar tendency with effective radius, 
supporting the robustness of our result. 

Dividing our sample into late and early-types according to 
the concentration index, we found that the colour-gradient of late- 
type galaxies change more rapidly with the galactic size than early- 
types. The corresponding slopes for samples SI 9, S20 and 525 are: 
-0.028, -0.025 and -0.024 for late-types and -0.004, -0.002 and 
-0.002 for early-types galaxies. The colour gradient of early-type 
galaxies, remains quite constant with galactic size. 

For the samples of galaxies under study, the colour gradient 
variation with R50 has been found to be independent of the amount 
of light used for obtaining it. Still, the deeming of colour gradi- 
ents could be partly due to calculating them with fewer annuli with 
surface brightness measurements. If this was the case, we would 
expect a correlation between the gradient error and R50(arcsec). 
However, no clear tendency is found, though the errors disper- 
sion clearly decreases for larger R50 (see Fig. lA2t . This tendency 
cannot explain alone the variation with R50 seen for colour gra- 
dients. However, this consideration brings us to point out the un- 
clear role that colour gradients could play in biasing global galactic 
properties obtained by extrapolating spectral c haracteristics from 
the SDSS galactic core s. In their studies, both [Scodeggio ( 2001) 
and lRoche et alj f2009) found out that the colour gradients should 
be taken into account for estimating the evolution of the colour- 
magnitude relation in early- type galaxies. 

Despite the large dispersion in colour gradients values, the 
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In the 520.5 sample, we find 128 galaxies that are within 
100 km/ s in redshift and 5 arcsec projected from another galaxy. 
A visual inspection of these galaxies reveals pairs of interacting 
galaxies, with some of them showing clear tidal tails. However, 
most of these close pairs appear just as red double nucleated struc- 
tures, similar to that shown in the inset of Fig. lA4l 

It is very interesting to see in Fig. lA4l that these galaxies tend 
to present bluer cores than what expected from the total sample. 
This result could point out that the star formation triggered by 
galactic interactions could be responsi ble for the observ ed bluer 
cores. Simulation such as the one from KobavashJ d2004h suggest 
that this phase will decay with time into a flatter colour gradient. In 
any case, these galaxies represent a percentage far below 1% of the 
initial 520.5 sample, and, therefore their contribution to increase 
the presence of blue nucleated galaxies is negligible at least in our 
sample. 
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Figure A4. Distribution of (g-i) colour gradient for all the galaxies within 
the S20.5 sample, solid line, and for only close pairs of galaxies within 
either 1 arcmin, dotted line, or 5 arcsec, dashed line, projected radius and 
100 km/ s in redshift. The inset shows an example of a pair of galaxies pair 
closer than 5 arcsec projected radius. Histograms are normalised to give 
unit area underneath. 




above results suggest that on average, galaxies with larger R50, 
measured in either arcsec or kpc, intrinsically tend to have steeper 
colour gradients. 



A3 Effect of having close pairs of galaxies 

Interacti ng galaxies can affect the colour gradients of each other 
(see e.g. lKo &Imll2005h . Moreover, colour gradient can be ill de- 
fined for close pairs of galaxies. Besides, the sky subtraction for 
close pairs of galaxies within the SDSS database may be problem- 
atic (John Lucey, private communication). Here we explore how 
the colour gradient distribute for close pairs of galaxies within our 
defined luminosity-threshold samples. 

Fig. lA4l show the distribution of colour gradients for galaxies 
that are within 100 km/ s in redshift and within a projected radius 
of either 1 arcmin or 5 arcsec. This constitutes a lower limit to the 
fraction of close galactic pairs, since this number is expected to 
increase when including the whole photometric SDSS sample^ 

We can observe that those close pairs with projected separa- 
tions of maximum 1 arcmin, present a distribution with a median 
value coincident with that of the total sample. These are only 5% of 
the total 520.5 sample. Within the studied redshift range, 1' corre- 
sponds to around 100 kpc, approximately half the typical Virial ra - 
dius of galaxies as massive as the Milky Way dCooper et alj|2009h . 
Thus, this will be an approximate limit to find interacting pairs of 
galaxies, which could affect each other colour gradients. Since no 
clear bias is found we do not have any further reason to consider 
that the colour gradients of close pairs of galaxies as far as 1 arcmin 
apart is ill defined. 



SDSS could not observe galaxies closer than 55 arcsec in the same plate 
and therefore close pairs are under represented in the spectroscopic sample. 
They were only observed in regions where spectroscopic plates overlapped. 



